arXiv:l507.06668v 1 [astro-ph.SR] 23 Jul 2015 


Astronomy & Astrophysics manuscript no. 26338arxiv ©ESO 2015 

July 27, 2015 


AMBER-NACO aperture-synthesis imaging of the half-obscured 
central star and the edge-on disk of the red giant L 2 Pup* * 

K. Ohnaka 1,2 , D. Schertl 2 , K.-H. Hofmann 2 , and G. Weigelt 2 


1 Universidad Catolica del Norte, Instituto de Astronomfa, Avenida Angamos 0610, Antofagasta, Chile 
e-mail: kl. ohnaka@gmail. com 

2 Max-Planck-Institut ftir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 
Received / Accepted 


ABSTRACT 

Aims. The red giant L 2 Pup started a dimming event in 1994, which is considered to be caused by the ejection of dust clouds. We 
present near-IR aperture-synthesis imaging of L 2 Pup achieved by combining data from VLT/NACO and the AMBER instrument of 
the Very Large Telescope Interferometer (VLTI). Our aim is to spatially resolve the innermost region of the circumstellar environment. 
Methods. We carried out speckle interferometric observations at 2.27 pm with VLT/NACO and long-baseline interferometric obser¬ 
vations with VLTI/AMBER at 2.2-2.35 pm with baselines of 15-81 m. We also extracted an 8.7 pm image from the mid-IR VLTI 
instrument MIDI. 

Results. The diffraction-limited image obtained by bispectrum speckle interferometry with NACO with a spatial resolution of 57 mas 
shows an elongated component. The aperture-synthesis imaging combining the NACO speckle data and AMBER data with a spatial 
resolution of 5.6x7.3 mas further resolves not only this elongated component, but also the central star. The reconstructed image reveals 
that the elongated component is a nearly edge-on disk with a size of ~180 x 50 mas lying in the E-W direction, and furthermore, that 
the southern hemisphere of the central star is severely obscured by the equatorial dust lane of the disk. The angular size of the disk is 
consistent with the distance that the dust clouds that were ejected at the onset of the dimming event should have traveled by the time 
of our observations, if we assume that the dust clouds moved radially. This implies that the formation of the disk may be responsible 
for the dimming event. The 8.7 pm image with a spatial resolution of 220 mas extracted from the MIDI data taken in 2004 (seven 
years before the AMBER and NACO observations) shows an approximately spherical envelope without a signature of the disk. This 
suggests that the mass loss before the dimming event may have been spherical. 
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1. Introduction 

Slow, but intense mass loss at the asymptotic giant branch (AGB) 
leads to the formation of thick circumstellar envelopes. When 
AGB stars evolve further to protoplanetary nebulae (PPNe) and 
planetary nebulae (PNe), striking bipolar lobes emerge, often ac¬ 
companied by collimated jets (Sahai & Trauger H998l) . Binarity 
is currently considered to be the most promising mechanism 
to shape PNe (e.g., De Marco 2009 ). but it is not yet under¬ 
stood how and at which evolutionary phase the bipolar structure 
emerges. High-resolution imaging of a few dusty AGB stars re¬ 
veals bipolar or more complex, clumpy structures (e.g., Weigelt 
et al. 119981 ; Monnier et al. 120001 Hofmann et al. 120011) . Some 
AGB stars show spiral structures, which strongly sugge sts the 
influence of a companion (e.g., Mauron & Huggins [20061 Maer- 
cker et al. 120121 Kim e t al. 120131 Ramstedt et al. 120141 Mayer 
et al. 12014 1; Decin et al. 2015 ). These observations suggest that 
the seed of the morphological change already exists in the AGB 
phase. 
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European Southern Observatory. Program ID: 074.D-0075(A), 074.D- 
0101(A), 074.D-0198(B), 088.D-0150(A/B), and 288.D-5041(A) 


The M5 giant L 2 Pup is a bright, nearby AGB star with a 
mass-loss rate of ~3 x 10 -7 M© yr _1 (Jura et al. 2002 ) that shows 
episodic, asymmetric dust formation. Magalhaes et al. ( 1986 ) in¬ 
terpreted temporally variable polarization in the optical as due to 
the growth and dissipation of dust grains in an asymmetric dust 
cloud. The 11.7 and 17.9 pm images of Jura et al. (12002 ) showed 
an asymmetric circumstellar dust envelope. Furthermore, Bed¬ 
ding et al. (2 0021 ) have detected a dimming event in the optical 
light curve of L 2 Pup starting in 1994, which suggests the ejec¬ 
tion of dust clouds. The brightness profiles reconstructed in the 
visible also lend support to the hypothesis that there are asym¬ 
metric dust clouds close to the star (Ireland et al. 2004 ). 

The circumstellar environment of L 2 Pup close to the star has 
recently been investigated in more detail. Kervella et al. ( 2014 , 
hereafter K14) observed L 2 Pup using the adaptive optics instru¬ 
ment NACO at VLT. Their images taken from 1.04 to 4.05 pm 
with a spatial resolution from 26 to 100 mas reveal a nearly edge- 
on disk lying in the E-W direction. The authors also showed 
that the observed images can be explained well by a flared disk 
model. The 2.30 and 3.74 pm images obtained by Lykou et 
al. ( 20151 . hereafter LI5) using the aperture-masking technique 
with NACO also show an elongated component. More recently, 
Kervella et al. (120151) have revealed a clear bipolar structure us¬ 
ing the VLT/SPHERE-ZIMPOL instrument. 
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Fig. 1 . uv coverage of our AMBER and NACO speckle observations of 
L 2 Pup. 


In this paper, we present near-IR interferometric obser¬ 
vations of L 2 Pup with an even higher spatial resolution of 
5.6x7.3 mas, combining VLT/NACO and the near-IR interfer¬ 
ometric instrument AMBER at the Very Large Telescope Inter¬ 
ferometer (VLTI). In addition, we also present a mid-IR image 
at 8.7 jum extracted from the mid-IR VLTI instrument MIDI. 


2. Observations 

The near-IR interferometric instrument VLTI/AMBER com¬ 
bines three telescopes and allows us to observe objects with a 
spatial resolution of ~3 mas (at 2 jun 1 ) with the currently avail¬ 
able baselines and a spectral resolution of 35, 1500, and 12000 
(Petrov et al. 2007 ). Our AMBER observations of L 2 Pup re¬ 
ported in this paper took place on 2011 December 23, 29, and 
30 (UTC) using two different telescope configurations (Al-Cl- 
D0 and D0-H0-I1), which provided projected baseline lengths 
from 15 to 81 m (Program ID: 088.D-0150A/B). The uv cover¬ 
age is shown in Fig. |T] We used the medium spectral resolution 
of 1500 in the wavelength region between 2.2 and 2.35 ju m, cov¬ 
ering the continuum and the CO first overtone bands. The Detec¬ 
tor Integration Time (DIT) was set to 100 ms in all AMBER data 
sets. We observed Sirius (a CMa, AIV, uniform-disk diameter = 
5.9 mas; Davis et al. 201 1) and Rigel (j8 Ori, B8Iae, uniform-disk 
diameter = 2.43 mas; Richichi et al. 120051) as interferometric cal¬ 
ibrators. 

Additionally, we carried out speckle interferometric observa¬ 
tions with VLT/NACO in Director’s Discretionary Time to ob¬ 
tain interferometric data at short baselines (Program ID: 288.D- 
5041A). This is crucial for the aperture-synthesis imaging of ob¬ 
jects with a very extended component. The NACO speckle ob¬ 
servations reported here occurred on 2012 March 21 with the 
IB2.27 filter, which is centered at 2.27 jum with a FWHM of 
0.06 jum. We recorded 398 frames with DIT = 109 ms using the 
S13 camera with a pixel scale of 13.2 mas. The window size was 
512 x 514 pixels, which resulted in a field of view of 678 x 678. 
We observed Canopus (a Car, A9II, angular diameter = 6.9 mas; 
Domiciano de Souza et al. 2008 1 as a calibrator and took 300 
frames with the same DIT as L 2 Pup. 

As complementary data, we also extracted mid-IR images 
at 8.7 jum (filter FWHM = 1.75 jum) from the acquisition im¬ 
age data obtained with the VLTI/MIDI instrument (Leinert et al. 
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2003 ). These images were taken to adjust the position of the tar¬ 
get on the detector. Thanks to the adaptive optics system for the 
MIDI operation with the 8.2 m Unit Telescopes (UTs) and the 
short DIT of 4 ms, the MIDI acquisition images are diffraction 
limited. In the ESO archive, there are MIDI data of L 2 Pup taken 
on 2004 December 31, 2005 January 3, 2005 February 28, and 
2005 March 3 (the interferometric data obtained on 2005 January 
3 (UTC) are presented in K14). The data taken on 2004 Decem¬ 
ber 31 with UT3 and UT4 (Program ID: 074.D-0198B) provide 
the best images of L 2 Pup and the calibrators, which are used as 
references of the point spread function (PSF), and we therefore 
present only these images in this paper. We reduced the data of 
s Cru and S Ori obtained on the same night to ascertain the PSF 
(Program IDs: 074.D-0101A and 074.D-0075A). As we show 
below (Fig. [4]), the quality of the PSF obtained from the data of 
s Cru is noticeably poorer than the image quality of L 2 Pup be¬ 
cause the former star is much fainter (Fu jjl m = 32.4 Jy) than the 
latter (F^jira = 2415 Jy). As a second PSF, we therefore used 
the Mira star S Ori because its mid-IR flux {F i 2 = 151 Jy) is 
much higher than £ Cru. Although S Ori was observed as a sci¬ 
ence target of interferometric observations, it is unresolved in the 
acquisition images at 8.7 yum, as we show below. The summary 
of our AMBER and NACO observations as well as the MIDI 
observations is given in Table. [I] 

We reduced the AMBER data with the amdlib ver. 3.0.3 
packag^] which is based on the P2VM algorithm (Tatulli et al. 
120071) . The product of the reduction is the visibility (or visibil¬ 
ity amplitude), closure phase, and differential phase. We binned 
the raw data in the spectral direction with a box car filter with a 
width of three pixels to increase the S/N ratio. The details of the 
reduction are described in Ohnaka et al. ( 2009 ). We excluded 
the data sets #2, #6-9, #13, and #15. They were taken under 
poor and variable seeing (^175) and/or short coherence time 
(^3.0 ms), which makes the absolute calibration of the visibil¬ 
ities unreliable. The NACO data were reduced with the bispec¬ 
trum speckle interferometry method (Weigelt fl9771 Lohmann et 
al. 119831 Hofmann & Wei gelt IT9861 . The result of the reduction 
of speckle data is the 2D visibility, the bispectrum on a huge 
number of baseline triangles within the maximum baseline of 
8 m (aperture of the telescope), and a reconstructed image of 
the object. For the combined image reconstruction with AMBER 
and NACO, we sampled visibilities and closure phases from the 
speckle data on randomly selected 855 baseline triangles within 
a maximum baseline length of 6 m (we avoided sampling be¬ 
tween 6 and 8 m because the data are noisy). 

For the reconstruction of an aperture-synthesis image from 
the combined AMBER and NACO speckle data, we used the 
image reconstruction package MiRA^] (Thiebaut 2008 ), which 
was also used for the image reconstruction of the red super¬ 
giants Betelgeuse and Antares in our previous works (Ohnaka 
et al. 120111 120131) . Given that the wavelengths covered by AM¬ 
BER and NACO overlap (2.2-2.35 and 2.27 jum, respectively), it 
can be justified to combine the AMBER and NACO data for the 
image reconstruction in the continuum. We used all wavelengths 
in the continuum region (2.2-2.29 /im) in the AMBER data to 
reconstruct one image because the object image is not expected 
to change noticeably in this narrow continuum region. Details of 
the image reconstruction procedure we adopted are described in 
Appendix [B] 

The reduction of the MIDI acquisition images is similar to 
the processing described in Ohnaka (120141 . After the sky sub- 


1 http: //w w w.j mmc .fr/data_proces sing_amber. htm 

2 http://cral.univ-lyon 1 .fr/labo/perso/eric.thiebaut/?Software/MiRA 
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Table 1 . Summary of the AMBER and NACO speckle observations and the archived MIDI data of L 2 Pup. The seeing, which was measured with 
a differential image motion monitor (DIMM), and the coherence time (r 0 ) were measured in the visible (see for details http://www.eso.org/gen- 
fac/pubs/astclim/paranal/seeing/adaptive-optics/ ). N f is the number of frames, and N exp is the number of exposures. 


# 

^obs 


PA 

Seeing 

To 

DIT x N f x N exp 

Remarks 


(UTC) 

(m) 

0 

O 

(ms) 

(ms) 




AMBER (A 1-Cl-DO): L 

2 Pup: 2011 Dec 23 (UTC) 


1 

01:48:07 

35.7/14.5/22.3 

-175/23/174 

0.9 

4.0 

100 x 500 x 5 


2 

02:33:57 

35.7/14.9/22.3 

-166/34/-180 

1.7 

2.5 

100 x 500 x 5 

Not used 

3 

03:11:28 

35.8/15.3/22.3 

-160/42/-174 

1.2 

3.0 

100 x 500 x 5 


4 

03:47:56 

35.8/15.6/22.2 

-153/50/-169 

0.9 

3.7 

100 x 500 x 5 


5 

04:24:20 

35.7/15.9/22.1 

-147/57/-164 

1.1 

3.5 

100 x 500 x 10 


6 

06:00:50 

34.4/15.8/21.4 

-132/74/-151 

1.7 

1.7 

100 x 500 x 10 

Not used 

7 

07:04:38 

32.5/15.1/20.5 

—123/85/—144 

2.2 

1.3 

100 x 500 x 15 

Not used 

8 

07:23:33 

31.8/14.7/20.2 

—121/88/—142 

2.2 

1.4 

100 x 500 x 10 

Not used 

9 

08:41:54 

27.4/12.8/18.2 

— 111/103/—134 

2.7 

1.1 

100 x 500 x 5 

Not used 

Cl 

01:19:52 

30.2/10.3/21.5 

-174/33/173 

0.8 

4.7 

100 x 1000 x 5 

Sirius 

C2 

02:16:06 

31.0/12.0/21.5 

-162/47.5/-178 

0.9 

4.1 

100 x 500 x 5 

Sirius 

C3 

03:30:29 

33.1/14.3/21.7 

-149/60/-167 

0.9 

4.1 

100 x 500 x 5 

Sirius 

C4 

04:06:45 

34.2/15.1/22.0 

-144/64/-162 

0.8 

4.7 

100 x 500 x 5 

Sirius 

C5 

04:50:11 

35.2/15.8/22.3 

-138/69/-157 

1.3 

2.7 

100 x 500 x 10 

Sirius 

AMBER (D0-H0-I1): 2011 Dec 29 (UTC) 

10 

06:46:01 

81.0/40.0/60.0 

114/—21/86 

0.9 

7.0 

100 x 500 x 8 


11 

07:38:25 

78.1/40.4/55.6 

124/-13/95 

0.9 

6.0 

100 x 500 x 5 


12 

08:20:26 

75.1/40.5/51.0 

134/—67/103 

0.9 

6.5 

100 x 500 x 5 


13 

08:56:34 

72.3/40.5/46.6 

142/-104/111 

1.0-1.7 

4.4 

100 x 500 x 5 

Not used 

C6 

06:23:18 

76.7/39.2/61.1 

107/-205/77 

1.0 

5.8 

100 x 500 x 5 

Sirius 

Cl 

07:10:38 

69.8/38.5/56.3 

113/-137/79 

2.0 

2.8 

100 x 500 x 1 

Sirius 

C8 

07:20:40 

73.1/40.6/49.0 

136/—44/105 

1.2 

4.6 

100 x 500 x 5 

Sirius 

AMBER (D0-H0-I1): 2011 Dec 30 (UTC) 

14 

02:28:00 

68.4/32.2/60.5 

67/-51/39 

0.7 

6.0 

100 x 500 x 5 


15 

05:23:07 

82.2/39.0/63.6 

100/-32/72 

1.5 

2.7 

100 x 500 x 5 

Not used 

C9 

00:09:57 

58.5/39.5/42.5 

95/-39/52 

1.2 

3.8 

100 x 500 x 5 

Rigel 

CIO 

01:10:28 

50.7/36.8/43.2 

83/-40/38 

1.1 

4.1 

100 x 500 x 5 

Sirius 

Cll 

01:52:41 

61.3/38.5/48.6 

87/-40/48 

0.9 

4.7 

100 x 500 x 5 

Sirius 

C12 

03:47:19 

81.4/39.0/64.0 

101/—29/73 

0.6 

7.6 

100 x 500 x 5 

Rigel 

NACO: 2012 Mar 21 (UTC) 

16 

02:30:14 

— 

— 

1.6 

3.0 

109 x 398 x 1 


C13 

02:38:11 

— 

— 

2.0 

3.0 

109 x 300 x 1 

Canopus 

MIDI (UT3-UT4): 2004 Dec 31 (UTC) 

17 

03:19:33 

— 

— 

1.2 

2.3 

4 x 1000 x 2 

Acquisition images 

18 

05:02:14 

— 

— 

1.7 

1.5 

4 x 1000 x 2 

Acquisition images 

C14 

01:47:07 

— 

— 

1.1 

2.6 

4 x 1000 x 2 

SOri 

C15 

02:37:28 

— 

— 

0.8 

3.6 

4 x 1000 x 1 

SOri 

C16 

07:09:13 

— 

— 

0.9 

2.7 

4 x 1000 x 2 

s Cru 


traction, we recentered each frame and added the frames. The 
sky-subtracted images of L 2 Pup show detector artifacts that are 
presumably caused by the high brightness of L 2 Pup. After re¬ 
moving these artifacts as described in Appendix [A| we added 
all frames in each data set of L 2 Pup. The calibrators are much 
fainter than L 2 Pup, and the residual of the sky subtraction is sig¬ 
nificant in some frames. Therefore, we only added the frames in 
which the residual of the sky subtraction is small, to obtain as 
clean a PSF as possible. The resulting image from each data set 
was derotated so that North is up and East is to the left, using 
the field rotation angle (Mathar [20061) . Finally, the derotated im¬ 
ages from all data sets were added. The quality of the images 
taken with UT3 is significantly lower than that taken with UT4, 
therefore we present the final image obtained from the UT4 data. 
We carried out neither the PSF subtraction nor the deconvolution 
because the images of L 2 Pup are saturated at the central peak. 


3. Results 

3.1. Near-IR AMBER-NACO aperture-synthesis image 

Figure [2] shows the images reconstructed from the NACO 
speckle data alone and from the combined NACO speckle and 
AMBER data. Both images are convolved with the following 
Gaussian beam. The beam size (FWHM) of the image recon¬ 
structed from the speckle data alone is 57 mas, which corre¬ 
sponds to the diffraction-limit of VLT at 2.27 jum. The beam 
size of the image reconstructed from the combined data is 
5.6 x 7.3 mas, which was derived by fitting the central peak of 
the dirty beam with a 2D elliptical Gaussian. Figure [3] shows a 
comparison of the visibilities and closure phases from the recon¬ 
structed image with the observed data. The reduced x 2 of the 
fit is 0.97. The image reconstructed from the speckle data alone 
shows a central star and an extended component elongated ap¬ 
proximately in the E-W direction with a size of ~300 x 200 mas. 
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Fig. 2. a: Image of L 2 Pup reconstructed from the NACO speckle data alone, b: Image reconstructed from the combined AMBER and NACO 
speckle data. 
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Fig. 3. Comparison of the observed interferometric data and the corresponding observables from the reconstructed image shown in Fig. [2^ (a 
color version of the figure is available in the electronic edition). The observed data are plotted with the red (black in the printed edition) symbols 
with the error bars shown in green (light gray), while the model is represented by the blue (dark gray) symbols, a: Visibilities, b: Enlarged view of 
the visibilities obtained by speckle interferometry with NACO. c and d: Enlarged views of the visibilities obtained with AMBER in two different 
spatial frequency ranges, e: Closure phases. 


The size of the extended component in the deconvolved images 
obtained by K14 and L15 in the K band is -200 x 100 mas, 
which roughly agrees with the measured size in our speckle im¬ 
age, which is convolved with the aforementioned beam. 

The image reconstructed from the combined NACO speckle 
and AMBER data, shown in Fig. [2 }), further resolves not only 
the elongated component seen in the speckle image, but also the 
central star. In this image, the extended component has a size of 
-180 x 50 mas. We checked that the NACO speckle + AMBER 
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image matches the NACO speckle-only image when convolved 
with the corresponding beam. Moreover, the figure shows that 
the southern half of the central star is severely obscured. This 
AMBER-NACO aperture-synthesis image directly shows both 
the edge-on disk and the half-obscured central star due to the 
equatorial dust lane of the disk. At the distance of 64 ± 4 pc 
(from the Hipparcos parallax of 15.61 ± 0.99 mas, van Leeuwen 
12007k the imaged size of the disk corresponds to (11.5 ± 0.8) x 
(2.8 ± 0.2) AU. 
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Fig. 4. 8.7 yt/m images of L 2 Pup and two PSF 
reference stars extracted from the MIDI acqui¬ 
sition data. In all panels, North is up and East is 
to the left, a: F 2 Pup with a minimum intensity 
of 0.1% of the central peak, b: L 2 Pup with a 
minimum intensity of 2% of the central peak to 
compare with the 11.7 jim image obtained by 
Jura et al. ( 20021 ). c: PSF reference star s Cru. 
The large arcs near the top and the bottom of the 
image are the residual of the background sub¬ 
traction. d: PSF reference star S Ori. 


Our combined aperture-synthesis image is qualitatively con¬ 
sistent with the disk model of K14. They also found out that the 
2-2.4 yum interferometric data of L 2 Pup taken in 2001 by the 
VLTI/VINCI instrument are better fitted with an ellipse with an 
axis ratio of 1.5 with the major axis at a position angle 106° (i.e., 
nearly in the E-W direction). This is naturally explained by the 
obscuration of the southern half of the central star by the edge- 
on disk, which makes the star appear to be elongated in the E-W 
direction. 

Our AMBER data cover the wavelength region of the CO 
first overtone bands. However, the speckle data were taken only 
in the continuum at 2.27 yum, not in the CO bands. This makes 
the image reconstruction more difficult, if not impossible, in the 
CO bands. Therefore, the analysis of the AMBER data in the CO 
bands will be presented in a separate, forthcoming paper. 

We note that the obscuration of the central star by a dust 
disk was also imaged for the eclipsing binary s Aur in the H 
band (1.50-1.74 /im) by Kloppenborg et al. ( 2010 1. In contrast to 
L 2 Pup, however, the disk of s Aur is located around an orbiting 
companion, and the disk itself is not visible at least in the H 
band. 

3.2. Mid-IR MIDI acquisition image 

Figure [4] shows the 8.7 /im images of L 2 Pup and the PSF refer¬ 
ence stars extracted from the MIDI acquisition data. The image 
of the first PSF reference s Cru (Fig. [4]:) shows noticeable resid¬ 
ual of the background subtraction, particularly near the edge of 
the field of view (large arcs near the top and bottom of the im¬ 
age). The image of the second PSF reference star S Ori is much 
better, and as Fig. [5] shows, the azimuthally averaged intensity 
profiles of two PSF reference stars agree, given the poorer qual¬ 


ity of the s Cru image. The FWHM of the intensity profiles of 
the PSF references indicates a a spatial resolution of 220 mas, 
which corresponds to the diffraction limit at 8.7 pm. The image 
of F 2 Pup (Fig. [4^) as well as the azimuthally averaged intensity 
profile (Fig. [5]) clearly shows an envelope much more extended 
than the PSF references. The envelope extends to an angular ra¬ 
dius of ~175 at the 0.1% intensity of the (saturated) central peak. 
Figure |4 ]d shows the image of F 2 Pup with an intensity cutoff of 
2% of the central peak, in the same manner as the 11.7 pm image 
presented by Jura et al. (2 0021 see their Fig. 2). Their 11.7 pm 
image taken February 2001 (i.e., about four years before the 
MIDI data) shows an envelope with a radius of ~1", while our 
MIDI acquisition image at 8.7 pm is much more compact with a 
radius of ~0'.'5 when shown with the same intensity cutoff This 
is probably because we probe the inner region of the envelope at 
8.7 pm than at 11.7 pm. For this reason, we cannot draw a con¬ 
clusion about time variations in the envelope between February 
2001 and December 2004. In marked contrast to the near-IR im¬ 
age presented above, the mid-IR image appears approximately 
spheric al w ithout a signature of the disk. We discuss this point 
in Sect. 14.31 

4. Discussion 

4.1. Comparison with the radiative transfer model of K14 

To examine whether our aperture-synthesis image of F 2 Pup 
is quantitatively consistent with the 2D radiative transfer (RT) 
disk model of K14, we compared the measured interferometric 

3 The central peak of the F 2 Pup image is saturated. If it were not for 
the saturation, the image with the 2% intensity cutoff would appear even 
more compact than shown in Fig. Ef 
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Fig. 5. Azimuthally averaged intensity profile of L 2 Pup (red solid line) 
and the PSF references s Cru (green dashed line) and S Ori (blue dotted 
line). The central region of L 2 Pup within O'.'l is saturated. 


observables—visibilities and closure phases—and the observed 
image with those predicted by their disk model. P. Kervella and 
M. Montarges kindly provided the best-fit model of their paper. 
From their model at 2.17 gm, which is the closest wavelength 
to our observations, we computed the visibilities and closure 
phases at the uv points of our NACO speckle and AMBER data. 

Figure [6] shows a comparison between the observed data 
and the RT model of K14. While the model reproduces the 
overall trend of the observed visibilities as a function of spa¬ 
tial frequency, the fit to the data is poor with the reduced x 2 = 
78.6. The model image at 2.17 gm convolved with the beam of 
our NACO+AMBER aperture-synthesis imaging (Fig.[6]f) shows 
that the disk appears fainter and more inclined (i.e., closer to 
edge-on) compared to the observed image shown in Fig. [2}) (the 
color scale is the same for the model and observed images). The 
intensity of the disk in our aperture-synthesis image is -7% of 
the central star, while it is -2% in the convolved image of the RT 
model. This means that the flux contribution of the star is higher 
in the RT model than in the observed image, which also explains 
that the model visibilities are higher than the observed data at 
spatial frequencies of (7 - 20) x 10 6 rad -1 . 

It is possible that the flux contribution of the central star was 
higher at the time of the observations of K14 than for our ob¬ 
servations. To examine this possibility, we changed the intensity 
of the star in the RT model image of K14. As Fig. [7] shows, the 
reduction of the stellar intensity by a factor of 0.6 can improve 
the match to the reduced^ 2 = 19.2. This result suggests that a 
time variation in the stellar flux by 0.55 magnitude at 2.2 pm 
between our AMBER and NACO observations (December 2011 
and March 2012) and the NACO observations of K14 (March 
2013) can at least partially explain the disagreement between 
the RT model of K14 and our data. The amplitude of the K -band 
light curve of F 2 Pup is ~0.5 magnitude (Whitelock et al. 12000b . 
Therefore, the variation in the stellar flux may be due to the pul¬ 
sation of the central star. However, there is still disagreement in 
the visibilities observed at the longest baselines (Fig. [7ji), and 
the disk appears too inclined compared to the reconstructed im¬ 
age. This implies time variations in the structure of the disk, for 
example, inclination angle and/or the flaring angle of the disk. 
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4.2. Binary companion 

Kervella et al. (2 0151 ) have recently found a companion at 
32.9 mas (2.1 AU) in west of the central star by adaptive op¬ 
tics polarimetric imaging in the visible with the VFT/SPHERE- 
ZIMPOF instrument. However, we do not see a signature of the 
companion in our NACO+AMBER image. A possible reason 
might be that the companion is too faint at 2.2 pm and/or our 
uv coverage is too sparse in the E-W direction, and therefore, 
we might have failed to detect the signature of the companion 
in the visibility or closure phase. Alternatively, the companion 
might have been at a different position at the time of our obser¬ 
vations. Kervella et al. d2015b estimated the orbital period to be 
1.4-4.6 years. It is possible that the companion was behind the 
near side of the inner rim of the disk and was invisible at the time 
of our AMBER and NACO observations. 

To set more quantitative constraints on the position and flux 
of the companion, we added a companion to the reconstructed 
image and examined the increase in the reduced^ 2 . The compan¬ 
ion was placed at random positions within 100 mas from the cen¬ 
tral star (corresponding to the AMBER field of view) and with 
a random flux contribution with respect to the star. We found 
out that the fit to the observed data is significantly poorer (re¬ 
duced x 2 > 2) for a companion with a flux contribution higher 
than 2.4% of the central star. This sets an upper limit on the flux 
of the companion at 2.2 pm. On the other hand, we can set no 
constraints on the position of the companion. 

The separation of the companion of 32.9 mas measured by 
Kervella et al. (12015b corresponds to 3.8 stellar radii, given the 
angular diameter of 17.5 mas derived by K14. This implies that 
the primary red giant star may be filling the Roche lobe and may 
possibly show deformation, as seen in the Roche-lobe-filling 
systems [5 Lyr and Algol = Per (Zhao et al. 120081 Baron et 
al. 2012 ). We examined this possibility by introducing elliptical 
deformation to the central star in the initial model of the im¬ 
age reconstruction. If the major-axis of the star (E-W direction) 
is longer than the minor-axis (N-S direction) by more than 40%, 
the reconstructed image cannot reproduce the observed data very 
well. The observed axis ratio of the Roche-lobe-filling star is 
-1.2 in the case of Lyr and Algol (Zhao et al. 2008 ; Baron et 
al. 2012 ). Therefore, our uv coverage is not sufficiently dense at 
long baselines to draw a definitive conclusion about the defor¬ 
mation. Aperture-synthesis imaging with a better uv coverage is 
necessary to probe the shape of the primary star. 

4.3. Dimming event and the disk formation 

As mentioned in Sect. [I] L 2 Pup started a dimming event around 
1994. If the dust clouds formed in 1994 move away at a radial 
velocity of -3 km s -1 (Winters et al. 2002 ), they should have 
traveled a distance of -11 AU ^170 mas by the time of our ob¬ 
servations. This is larger than the angular radius of the disk seen 
at 2.2 pm, but comparable to the radius of the disk seen at 3-4 pm 
(K14; LI5). This suggests that the disk may have been formed 
by the dust clouds that have been ejected since the onset of the 
dimming event. The dust clouds ejected at the very beginning 
are perhaps already too far away from the star and, therefore, 
too cold to be detected at 2.2 pm. We note, however, that this 
argument is based on the assumption that the dust clouds move 
away radially. The 3.74 pm images obtained by K14 and LI5, the 
4.05 pm image of K14, and the visible image of Kervella et al. 
( 2015 ) show a spiral structure, which is presumably caused by 
the companion. If the dust clouds move along the spiral struc¬ 
ture, the radial distance traveled by the dust clouds is smaller 
























Ohnaka et al.: AMBER-NACO imaging of the half-obscured central star and the edge-on disk of L 2 Pup 


a. 


a 

H0.5 h 
& 


• Observation 

• K14 model 
Xred = W.6 


.'I 

I 




10 ' 


2x10' 


3x10- 


Spatial frequency (rad ) 



Spatial frequency (rad 1 ) 


5x10° 10' 1.5x10 7 

Spatial frequency (rad -1 ) 


2x10' 



100 



-100 


100 


50 0 -50 

RA offset (mas) 


-100 


Fig. 6. Comparison of our NACO speckle and AMBER observations with the RT model of K14 (a color version of the figure is available in the 
electronic edition). In panels a-e, the observed data are plotted with the red (black in the printed edition) symbols with the error bars shown in 
green (light gray), while the model is represented by the blue (dark gray) symbols, a: Visibilities, b: Enlarged view of the visibilities obtained by 
speckle interferometry with NACO. c and d: Enlarged views of the visibilities obtained with AMBER in two different spatial frequency ranges, e: 
Closure phases, f: Model image in the same color scale as the reconstructed images in Fig. [2] 
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Fig. 7. Comparison of our NACO speckle and AMBER observations with the modified RT model of K14, shown in the same manner as in Fig. [6] 
A color version of the figure is available in the electronic edition. 


than estimated above. Therefore, kinematical information of the 
spiral structure is necessary to confirm the connection between 
the formation of the disk and the onset of the dimming event. 


well as the 11.7 n m and 17.9 gm images presented in Jura et al. 
( 2002) do not show a clear signature of the disk in this direction. 
Our MIDI acquisition image appears approximately symmetric. 
The mid-IR images of Jura et al. (120021) are characterized by 
an elongation at a position angle of -135° and a distinct blob 


While the disk seen in the near-IR is oriented roughly in the 
E-W direction, our mid-IR image at 8.7 gm shown in Fig. [4] as 
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at 225°. This makes the object appear slightly elongated in the 
E-W direction with a size of - 2'.'5 x 270, but the elongation 
is much less pronounced than seen in the near-IR. Their spatial 
resolution is 0'/47 and 0Y49 at 11.7 and 17.9 pm, respectively, 
which is significantly smaller than the angular size of the object, 
and therefore, the lower spatial resolution cannot fully account 
for the difference in the morphology. This means that the disk 
only exists in the innermost region of the circumstellar environ¬ 
ment, while the outer region is approximately spherical at least 
at the time of the observations of Jura et al. ( 2002 ) and the MIDI 
observations. 

The mid-IR imaging observations of Jura et al. (120021) took 
place in February 2001, while the MIDI data were obtained on 
December 2004. The dust clouds ejected at the onset of the dim¬ 
ming event (around 1994) should have traveled an angular dis¬ 
tance of -70 mas (disk diameter of -140 mas) and -100 mas 
(disk diameter of -200 mas) at the time of Jura et al.’s obser¬ 
vations and the MIDI observations, respectively. These angular 
displacements are smaller or just comparable to the angular res¬ 
olution of these mid-IR observations, and therefore may have 
been too small to detect. Furthermore, while we assumed that the 
dust clouds move away radially, they may actually move along 
the spiral structure, in which case the radial distance would be 
smaller. This might also have been a reason for the non-detection 
of the disk in the mid-IR observations. The circumstellar enve¬ 
lope imaged by these mid-IR observations most likely represents 
a nearly spherical mass loss that was already present before the 
formation of the disk. 

If we assume that the dust clouds move away radially, the 
dust clouds ejected at the beginning of the dimming event should 
have reached an angular distance of -200 mas by now (i.e., the 
disk diameter of -400 mas), which is resolvable with mid-IR in¬ 
struments of the current 8-10 m telescopes. If the dust clouds 
move along the spiral structure, the angular size of the enve¬ 
lope should be smaller. Therefore, new mid-IR imaging would 
be interesting to examine whether or not the morphology of the 
circumstellar envelope on larger spatial scales has changed due 
to the formation of the disk and whether the measured size is 
consistent with the radial motion of the dust clouds or not. 

5. Concluding remarks 

We have obtained the aperture-synthesis image of both the disk 
and the central star of F 2 Pup by combining the VFT/NACO 
speckle data and VFTI/AMBER long-baseline interferometric 
data. The obtained image shows that the disk has a size of -180 
x 50 mas and the southern half of the central star is severely 
obscured by the equatorial dust lane of the nearly edge-on disk. 
This agrees with the recent AO and aperture-masking images and 
modeling. The observed size of the disk is consistent with the an¬ 
gular distance that the dust clouds formed at the beginning of the 
dimming event (-1994) should have traveled by the time of our 
observations, if we assume that the dust clouds move away ra¬ 
dially. This lends support to the hypothesis that the formation of 
the dust disk is responsible for the dimming event. However, the 
kinematical information of the spiral structure is important for 
establishing the connection between the formation of the disk 
and the dimming event. AEMA observations of molecular lines 
are ideal for this goal. 

The mid-IR image extracted from the MIDI data taken about 
seven years before our AMBER and NACO observations shows 
an approximately spherical envelope, which most likely results 
from the mass loss occurring before the disk formation. High- 
resolution mid-IR imaging is useful for studying possible mor¬ 
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phological changes in the circumstellar environment on larger 
spatial scales since the onset of the dimming event. 
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Appendix A: Removal of detector artifacts from 
MIDI acquisition images 

As Fig. | A. lfr shows, noticeable vertical stripes are present in the 
sky-subtracted frames of the MIDI acquisition data of L 2 Pup. 
These detector artifacts presumably result from the high bright¬ 
ness of L 2 Pup (the central region within O'.'l is saturated). No 
vertical stripes appear in the images of the PSF reference stars. 
The intensity of the vertical stripes is different in the upper and 
lower regions with respect to the star, and it therefore is neces¬ 
sary to remove the stripes in the upper and lower regions sepa¬ 
rately. In each region, we removed the stripes as follows. First, 
for each column affected by the vertical stripes, we computed the 
median of the pixel values from the rows sufficiently far away 
from the star. Then we subtracted this median from all pixels in 
the column. We carried out this procedure for all columns that 
were affected by the stripes. Figure |A.1| demonstrates that the 
image is nearly free from the stripes after this procedure. 
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Fig. A.l. Removal of the detector artifacts in 
the MIDI acquisition images of L 2 Pup. a: One 
of the sky-subtracted frame of L 2 Pup, showing 
vertical stripes in the columns near the center, 
b: Same frame after the removal of the vertical 
stripes. The color scale in the diffraction core is 
saturated to clearly show the stripes. 


Appendix B: Image reconstruction 

The observed visibilities plotted in Fig. [3^ show a steep drop 
at low spatial frequencies ^ 3 x 10 6 rad - ^^ baselines shorter 
than 6 m), which suggests a very extended component. The vis¬ 
ibilities observed at spatial frequencies from ~6 x 10 6 rad 1 to 
10 7 rad -1 (baselines from 13 to 22 m) may appear reminiscent 
of the first and the second visibility lobe expected from a uni¬ 
form disk or limb-darkened disk. However, the visibility from a 
uniform disk without an extended component is 0.13 (or lower 
for a limb-darkened disk) in the extrema of the second visibil¬ 
ity lobe, much lower than the observed values of ~0.3. With an 
extended component as revealed by the speckle data, the visi¬ 
bility would be even lower. We first attempted to explain these 
observed data using geometrical models. While simple geomet¬ 
rical models may not fit the data completely, they are useful for 
characterizing the approximate geometry of the object and can 
also be used as an initial model for the image reconstruction. 


Appendix B. 1: Geometrical model: uniform disk + elliptical 
Gaussian 


We tried to fit the data with a geometrical model consisting 
of a uniform-disk-like central star and an elliptical Gaussian. 
The free parameters are the uniform-disk diameter of the cen¬ 
tral star (0*), the fractional flux contribution of the central star 
/*, the widths of the elliptical Gaussian along the major and mi¬ 
nor axes cr major and cr minor (the elliptical Gaussian is given by 


^-((A/cr ma j or ) 2 +0;/(T m i n or) 2 )^ anc [ p 0S iti 0 n angle of its major axis 
PA (measured from North to East). We searched for the best-fit 
model by varying 0* = 8 ... 22 (mas) with A0* = 2 (mas), /* = 
0.1 ... 0.7 with A/* = 0.05, <x ma j 0 r = 30 ... 100 (mas) with Acr ] 


' major 


= 10 (mas), o ~minor = 10 ... 50 (mas) with Acrminor = 10 (mas), and 
PA = 70°... 100° with APA = 5° (the PA was limited to this range 
from the elongation of the image reconstructed from the speckle 
data alone). The best-fit model, which is plotted in Fig. ED is 
characterized by 0* = 12 mas, /* = 0.3, <x ma j or = 50 mas, cr m i nor 


20 mas, and PA = 95° with a reduced^ 2 of 27.2. Figure 


B.l 


reveals that the best-fit model cannot reproduce the observed vis¬ 
ibilities at spatial frequencies of (0.6 - 2.0) x 10 7 rad -1 (baseline 
length = 13-44 m). Most of the observed visibilities at these 
spatial frequencies are noticeably higher than predicted by the 
model. This means that there is some sharp structure that is not 
seen in the smooth Gaussian model. 


Appendix B.2: Geometrical model: uniform disk + elliptical 
ring 


We then tried a model consisting of the uniform-disk-like central 
star and an elliptical ring, because a ring gives rise to visibilities 
much higher than a uniform disk or limb-darkened disk at long 
baselines. The free parameters are the uniform-disk diameter of 
the central star 0*, the fractional flux contribution of the central 
star /*, the semi-major and semi-minor axes of the elliptical ring 
CRmajor and RmmoT, respectively), and its position angle PA. The 
width of the ring was set to be 10% of its radius at each position 
angle. We searched for the best-fit model by varying 0* = 8 ... 
22 (mas) with A 0* = 2 (mas), /* = 0.1 ... 0.6 with A/* = 0.05, 
^major — 50 ... 100 (mas) with A7? ma j 0 r — 10 (mas), R m'm nr — 10 
... 50 (mas) with A 7? ma j o r = 10 (mas), and PA = 70° ... 100° with 
APA = 5°. Figure B.2 shows a comparison of the best-fit ring 
model with the observed data. This model is characterized by 
0* = 16 mas, /* = 0.6, R ma j or = 60 mas, 7? m i nor = 40 mas, and 


PA = 95° with the reduced^ 2 = 54.8. While the visibilities at 
(0.6 - 1) x 10 7 rad -1 from this model is as high as or even higher 
than the observed data, the fit to the speckle data and the data at 
longer baselines is worse than the above star + Gaussian model 
as shown by the worse reduced^ 2 . 


Appendix B.3: Geometrical model: obscured uniform disk + 
elliptical Gaussian 

We found out that the observed data are much better reproduced 
if the southern half of the aforementioned uniform-disk-like cen¬ 
tral star + elliptical Gaussian model is obscured. In this model, 
the intensity of the uniform-disk star + elliptical Gaussian model 
I 0 (x,y) (x and y are the coordinates on the sky with the origin at 
the central star) is modified as follows: 

I(x,y) = / 0 (x,y)x 

where s is a parameter to smoothen the obscuration edge. We set 
s to be 0.2, which decreases the intensity to zero over ~2 mas 
in the y direction. The free parameters are the uniform-disk di¬ 
ameter of the central star, the fractional flux contribution of the 
central star /*, the widths of the elliptical Gaussian along the 
major and minor axes <x ma j or and cr m i n0 r, and the position angle 
of its major axis PA. We searched for the best-fit model by vary¬ 
ing 0* = 8 ... 22 (mas) with A0* = 2 (mas), /* = 0.2 ... 0.6 with 
A/* = 0.05, <r ma jor = 30 ... 100 (mas) with A<x ma j or = 10 (mas), 
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Fig. B.l. Best-fit model consisting of a uniform-disk central star and an elliptical Gaussian shown in the same manner as Fig. [6] 
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Fig. B.2. Best-fit model consisting of a uniform-disk central star and an elliptical ring shown in the same manner as Fig. [6] 
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cr minor = 10 ... 50 (mas) with A<x m inor = 10 (mas), and PA = 70° 
... 100° with APA = 5°. 

The best-fit model is characterized by 0* = 20 mas, /* = 
0.45 , or m ajor = 70 mas, cr minor = 30 mas, and PA = 85°. As 
Fig. |B.3| shows, this model can reproduce the observed visibil¬ 
ities much better than the above two models. The uncertainties 
in /*, cr major? t^minor, and PA are ±0.05, ±10 mas, ±10 mas, and 
±5 °. The uniform-disk diameter of the central star is in the range 
between 18 and 20 mas. This agrees with the 17.5 mas derived 
by K14, given the difference in the data and the model used by 
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them and us. The reduced x 2 of this model is 16.2, which is 
much better than the star + Gaussian or star + ring models, but 
still much larger than 1. This is because the fit to the visibilities 
observed at the longest baselines (spatial frequencies higher than 
~2 x 10 7 rad -1 ) and to the observed closure phases are not sat¬ 
isfactory. However, this disagreement can be due to small-scale 
structures that are not included in our geometrical model, but can 
be modeled by the image reconstruction. 





















Ohnaka et al.: AMBER-NACO imaging of the half-obscured central star and the edge-on disk of L 2 Pup 




100 



-100 


100 


50 0 -50 

RA offset (mas) 


-100 


Fig. B.3. Best-fit model consisting of a half-obscured uniform-disk central star and an elliptical Gaussian shown in the same manner as Fig. [6] 


Appendix BA: MiRA parameters 

We used the best-fit half-obscured star + elliptical Gaussian 
model as the initial model. The regularization scheme of the 
maximum entropy method was adopted, with the prior being the 
Gaussian with the same widths as the initial model. The degree 
of regularization was set to p = 10 3 (see Thiebaut 1 200811 for de¬ 
tails of the regularization scheme and the definition of p). We 
also reconstructed the image with the total variation regulariza¬ 
tion scheme, but the image shows no noticeable differences. 
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